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The QCD expectations concerning the deep inelastic lepton - hadron scattering at low x and 
their phenomenological implications for HERA are summarised. Theoretical predictions for 
the structure function F2(x,Q ) based on the leading logl/x resummation are presented and 
compared with the results obtained from the Altarelli-Parisi equations. Theoretical predictions 
are confronted with the recent data from HERA. The role of studying the hnal states in deep 
inelastic scattering for revealing the dynamics at low x is emphasised and some dedicated 
measurements like deep inelastic plus jet events, transverse energy flow and dijet production 
in deep inelastic scattering are discussed. 



In this talk we shall give an overview 
of the QCD expectations concerning the 
small x behaviour of deep inelastic lepton- 
hadron scattering and will discuss their 
phenomenological implications for HERA. 
Besides discussing the deep inelastic struc- 
ture functions we shall also consider spe- 
cific measurements like deep inelastic scat- 
tering accompanied by the energetic jet, 
transverse energy flow and dijet produc- 
tion in deep inelastic scattering which 
should test the QCD predictions at small 
x more directly and unambiguosly. 

The dominant role at small x is played 
by the gluons and so at first we consider 
the small x limit of the gluon distribu- 
tion. The gluons are not, of course, di- 
rectly probed in deep inelastic lepton scat- 
tering, yet they contribute indirectly to 
this process through the g — > qq transi- 
tions. 



The relevant framework for calculating 
the parton distributions in perturbative 
QCD in the small x limit is the leading 
log(l/a;) (LLl/x) approximation. It cor- 
responds to the sum of those terms which 
contain the maximal power of ln(l/x) 
at each order of the perturbative expan- 
sion. The basic quantity in this approxi- 
mation is the unintegrated gluon distribu- 
tion f(x,k 2 ) which is related in the follow- 
ing way to the conventional (scale depen- 
dent) gluon distribution g(x,Q 2 ): 

rQ 2 dk 2 

xg(x,Q 2 )=j o -^f^k 2 ) (1) 

In the LLl/x approximation the uninte- 
grated gluon distribution satisfies the fol- 
lowing equation |3]-||: 
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which is called in the literature the 
Balitzkij-Fadin-Kuraev-Lipatov (BFKL) 
equation. 

It corresponds to the sum of ladder 
diagrams but unlike the leading logQ 2 
approximation which gives the Altarelli- 
Parisi evolution equations || the trans- 
verse momenta of the gluons are not or- 
dered along the chain. The kernel Kl 
of the eq.((2|) contains both the real gluon 
emission term as well as the virtual cor- 
rections. The former corresponds to the 
term proportional to f(x, k' 2 ) while the 
latter correspond to terms proportional to 
f(x,k 2 ) in the integrand on the rhs. of the 
eq.(||). The virtual corrections correspond 
to the " gluon reggeisation" Q f| (or to the 
"non-Sudakov" form factor ||). The 
variable(s) k 2 (k' 2 ) denote the transverse 
momenta squared of the gluons along the 
ladder. The parameter kg is the infrared 
cut-off which is necessary if the running 
coupling constant effects are taken into ac- 
count. More general treatment of the in- 
frared region in the BFKL equation than 
simply imposing the lower limit cut-off A; 2 
has been discussed in f7j-|§. 

When the running coupling effects are 
neglected (i.e. when one sets a s (k 2 ) = a s ) 
and when fcg = then the equation (|2|) 
can be solved analytically and the leading 
small x behaviour of its solution is given 
by the following formula: 

f(x,k 2 )~ 
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where the Riemann zeta function £(3) = 
1.202. The parameter k is of nonpertur- 
bative origin and is fixed by the boundary 
condition for f(x, k 2 ) at x = xq ||. 

The solution @ of the BFKL equa- 
tion follows from the solution of the 
corresponding equation for the moment 
f(n,k 2 ) 



f(n,k 2 ) 



dxx n ' 2 f{n,k' 2 ) (6) 



of the distribution f(x,k 2 
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where f°(n, k 2 ) is the moment of the suit- 
ably defined inhomogeneous term. This 
equation is diagonalised by the Mellin 
transform and its solution for the Mellin 
transform /(n, 7) of the moment f(n,k 2 ) 
is: 

/ (n, 7 ) 
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where 

K L {i) = 2*(1) - *( 7 ) - #(1 - 7) (9) 

is the Mellin transform of the kernel of 
the equation (^). The function ^(z) is 
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the logarithmic derivative of the Euler T 
function. The Mellin transform f(n, 7) is 
defined as below: 



/(>»> 7) 



dk 2 (k 
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f(n,k 2 ) (10) 



The poles of f(n, 7) in the 7 plane define 
the gluon anomalous dimensions 7(71, a s ). 
The leading twist anomalous dimension 
j LT (n,a s ) controls the k 2 dependence of 
f(n,k 2 ) at large k 2 



f(n,k z )~(k z ) 



2\~f LT (n,a s ) 



(11) 



It follows from the eq. (|8|) that the 
anomalous dimensions are the solutions of 
the following equation: 



3a s 



7r(n — 1) 



K L (j(n,a s )) = 1 
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The solution of this equation allows to ob- 
tain the gluon anomalous dimensions as a 
power series of 



3a s 



7r(n— 1) 



[llj, |12|. For the 



leading twist anomalous dimension this 
power series corresponds to the leading 
logl/z expansion of the splitting function 
P gg (z,a s ) which appears in the evolution 
equation for the gluon distribution. The 
leading twist quark anomalous dimension 
which includes the resummation of the 
powers of 7r ^f 1 ) has recently been dis- 
cussed in the ref. [13|. 



The exponent A controlling the small x 
behaviour of f(x, k 2 ) (cf. equations (||, |])) 
is 

A = ^K L {\/2) (13) 

7T 

The leading twist anomalous dimension 
has a branch point singularity at n = 1+A. 
We also have (cf. eqs. (12, PD): 



7 LT (n = 1 + \,a s 



(14) 



The following properties of the solution 
of the BFKL equation summarised in the 
formula (|3|) should be noted: 

(1) It exhibits the Regge type x~ A in- 
crease with decreasing x where the expo- 
nent A given by eq. (4) can have poten- 
tially large magnitude ~ 1/2 or so. The 
quantity 1 + A is equal to the intercept of 
the so called BFKL Pomeron which corre- 
sponds to the hard QCD interactions. 

Its potentially large magnitude (~ 1.5) 
should be contrasted with the intercept 
a soft ~ 1.08 of the (effective) "soft" 
Pomeron which has been determined from 
the phenomenological analysis of the high 
energy behaviour of hadronic and photo- 
production total cross-sections [14]. 

(2) It exhibits the (A; 2 ) 1 / 2 growth with 
increasing k 2 modulated by the Gaussian 
distribution in ln(k 2 ) with its width in- 
creasing as ln l / 2 (\/x) with decreasing x. 
The Gaussian factor reflects the diffusion 
pattern of the BFKL equation §, 
The increase of the function f(x,k 2 ) as 
{k 2 ) 1 / 2 is due to the fact that the leading 
twist anomalous dimension is equal to 1/2 
at the BFKL singularity (cf. eq. (0)). 
This shift of the anomalous dimension is 
the result of the (infinite) LL(l/x) re- 
summation and should be contrasted with 
the anomalous dimension calculated per- 
turbatively retaining only finite number of 
terms [|TJ]. 

If the running coupling effects are taken 
into account then the leading small x be- 
haviour is: 



f(x,k 2 ) 



~ X 



(15) 



The exponent A has to be now calculated 
numerically and is dependent upon the in- 
frared cut-off kl |], g, ||, |1|| . More recent 
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discussion of the BFKL equation is given 

MM- 



in 



The validity of the BFKL equation is, 
in principle, restricted to the region where 
a s ln(l/x) ~ 0(1) and a s ln{Q 2 /Ql) « 1 
where Qq is some moderately large scale 
(Qq >> A 2 ). Possible extension of the 
BFKL equation beyond this region is pro- 
vided by the Marchesini equation [|| which 
treats both (large) logarithms ln(l/x) and 
ln{Q 2 /Qq) on equal footing. In the region 
of large values of x this equation becomes 
equivalent to the Altarelli-Parisi evolution 
equations. 

For xg(x, Q 2 ) growing as x~ x the trans- 
verse area S(x, Q 2 ) occupied by the gluons 



S(x, Q 2 ) = const xg(x, 



a s (Q 2 ) 
Q 2 



(16) 



can, for sufficiently small value of x and for 
fixed Q 2 become comparable to the trans- 
verse area of a hadron Sh = kR 2 ^ where 
Rh is the hadronic radius. When this hap- 
pens (and in fact before this happens), the 
gluons can no longer be treated as free par- 
tons and their interaction leads to screen- 
ing (or shadowing) effects k| |19| . The 
main effect of shadowing is to tame the 
indefinite increase of parton distributions. 
One finds instead that at sufficiently small 
values of x and/or Q 2 the gluon distri- 
butions approach the so-called saturation 
limit xg sat {x,Q 2 ) |, |] 



xg sa t(x,Q' 



const 



R\Q 2 (17) 



In some models |20[j the saturation limit 
contains some remnant weak x depen- 
dence. The most dramatic effect is 
the linear scaling violation exhibited by 
xg S at(x,Q 2 )- 



If one assumes that the gluons are 
not distributed uniformly within a hadron 
but are concentrated around the "hot- 
spots" pi], 20 1 having their radius Rh. s . 
much smaller than the hadronic radius 
Rh then the shadowing effects are ex- 
pected to be stronger. The saturation 
limit xg sa t(x, Q 2 ) is then controlled by the 
radius Rh.s. and not by Rh- 

The shadowing effects modify the BFKL 
equation by the non-linear terms [|], [|, [0| 
2§: 



— x- 



df(x, k 2 
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The equation ( |18|) is called in the literature 
the Gribov, Levin, Ryskin (GLR) equa- 
tion. 

The second term in the right hand side 
of the eq.([l^) which is quadratic in the 
gluon distribution describes the shadow- 
ing effects. They correspond to the QCD 
diagrams where two gluonic ladders merge 
into one. The parameter R describes the 
size of the region within which the gluons 
are concentrated. 

The GLR equation has been derived as- 
suming that the gluons within a hadron 
are uncorelated and that the interaction 
between the gluonic ladders can be ne- 
glected. These assumptions have recently 
been shown to be unjustified p4| . The in- 
teraction between gluonic ladders can be 
approximmately accounted for by intro- 
ducing the appropriate enhancement fac- 
tor in the nonlinear term in the eq. ( Jl8| ) 

in- 

It turns out that in the region of x and 
Q 2 which may be relevant for HERA the 
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shadowing effects are still relatively weak 
at least for R ~ ^GeV^ 1 (T^j. The possi- 
bility of detecting the nonlinear shadowing 
terms in HERA has recently been discused 
in detail in the ref. [26]. 

In Fig.l we show the regions in the 
(lnl/x, ln(Q 2 /A 2 )) plane where various 
dynamical effects and approximations dis- 
cussed above are expected to play the 
dominant role. 

We shall now discuss quantitative 
predictions for the structure function 
F2(x, Q 2 ) at small x which will be directly 
based on the solution of the BFKL equa- 
tion ||, [27], ^8|. (For the first attempt to 
determine the structure function ^(x, Q 2 ) 
from the solution of the BFKL equation 
see the ref. ]p9fl ). The relevant diagrams 
are shown in Fig. 2 and their contribution 
to the structure function ^(x, Q 2 ) can be 
written in the following factorisable form: 



F 2 (x,Q 2 



dx f dkrp 



(19) 

where x/x' is the longitudinal momentum 
fraction carried by the gluon which cou- 
ples to the qq pair. Similar representation 
can be written for the longitudinal struc- 
ture function Fl(x,Q 2 ). The function F% 
which denotes the quark box contribution 
to the photon-gluon subprocess shown in 
Fig.l is given in ref. [27]. The function 
f(x/x',k^) is the unintegrated gluon dis- 
tribution corresponding to the sum of lad- 
der diagrams in the lower part of the di- 
agrams shown in Fig. 2. In the leading 
log(x'/x) approximation this distribution 
is given as the solution of the BFKL equa- 
tion (0). The (x/x')~ x behaviour of the 



unintegrated gluon distribution / (see the 
equations [l5])) generates the singular 
x~ x behaviour of the structure function 
F2(x,Q 2 ). We may also study the effects 
of shadowing in the structure functions by 
using the function / which is the solution 
of the GLR equation (||). 

The Q 2 dependence of the structure 
function F2(x,Q 2 ) reflects the k\ depen- 
dence of the function f(x,k^) p8fl . Thus, 
for instance, the increase of the uninte- 
grated gluon distribution as (fc^) 1//2 with 
increasing k\ (see the eq. (||)) implies in- 
crease of the structure function F2(x,Q 2 ) 
as {Q 2 ) 1 / 2 with increasing Q 2 . 

The formula ((l^) is an example of the 
kx factorisation theorem which is the basic 
tool for calculating the observable quanti- 
ties (i.e. in this case the structure func- 



tions) at small x [30 1 . Its connection with 
the more conventional collinear factorisa- 
tion has recently been discussed in detail 
in the ref. [31]. 



In Fig. 3 we summarise the predictions 
for the structure function F 2 {x, Q 2 ) based 
on the formula ([n]) [||, |27]] (the curves 
marked as AKMS) and confront them with 
the recent data from HERA fl3q, 3TL |38l. 



We also compare these predictions with 
structure functions which were obtained 
within the NLO Altarelli-Parisi formalism. 
Let us recall that in the latter case the 
small x behaviour depends upon the phe- 
nomenological input distributions at the 
reference scale and is therefore not con- 
strained by a theory. One can mimick the 
BFKL behaviour at small x by assuming 
the x~ x type extrapolation of the input 
parton distributions as it has been done in 
the refs. Q ||] (the curves MRS (A), 
f34}] ). In the next-to- leading approxima- 
tion one neglects, of course, the effects 
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of the leading logl/x resummation in the 
splitting (and coefficient) functions which 
are automatically taken care of in the cal- 
culation based on the fey factorisation with 
the function / obtained from the solution 
of the BFKL equation. Systematic study 
of those effects within the Altarelli -Parisi 
formalism has however shown that their 
role is relatively small provided the start- 
ing distributions are sufficiently singular 
at small x [|l2|. This result explains possi- 
ble origin of similarity between the AKMS 
and MRS curves. 

In the plots in Fig. 3 which correspond 
to Q 2 = 15 GeV 2 and to Q 2 = 30 GeV 2 
we also show predictions obtained from 
the "dynamical model" of parton distri- 



butions [35]. In this model the parton dis- 
tributions are generated radiatively from 
the valence-like input at the very low ref- 
erence scale Ql ~ 0.25GeV 2 . The strong 
increase of F2(x, Q 2 ) with decreasing x at 
the HERA range (i.e. for Q 2 ~ WGeV 2 ) 
comes now from the relatively large evolu- 
tion length i{Q 2 ,Ql) 



Q 2 dq 2 3a s (q 2 



TT 



(20) 



IQI 1* 

since in this case F2{x,Q 2 ) behaves at 
small x as below: 



F 2 (x,Q 2 ) ~ exp \2^(Q 2 ,Q 2 )ln(l/x) 

(21) 

The structure function Fzix, Q 2 ) which 
is generated radiatively in the lead- 
ing logQ 2 approximation from the non- 
singular input distributions at small x, 
exhibits scaling in the two variables 
y/ln(l/x)lnlnQ' 2 and y / /n(l / 'x) /InlnQ' 2 . 
Exploration of this scaling property has 
recently been advocated possible dis- 
criminator between the BFKL-motivated 



or purely radiatively generated structure 
functions and evidence has been found for 
supporting the latter p9[ . Above scaling 
property of structure functions may how- 
ever be modified if the leading log(l/x) 
terms are included in the splitting and co- 
efficient functions. These terms play sig- 
nificant role in the case of the non-singular 



starting parton distributions [12|. 

We notice from Fig. 3 that in the region 
of x and Q 2 relevant for HERA measure- 
ments the structure function F2(x,Q 2 ) 
which was based on the solution of the 
BFKL equation is similar to that which 
was obtained within the Altarelli-Parisi 
evolution equations formalism in the NLO 
approximation. The inclusive quantity 
like F 2 is not therefore the best discrim- 
inator for revealing the dynamical details 
at low x. 

One may however hope that this can be 
provided by studying the structure of the 
final states at small x. It is expected in 
particular that absence of transverse mo- 
mentum ordering which is the characteris- 
tic feature of the BFKL dynamics should 
reflect itself in the less inclusive quantities 
than the structure function F 2 (x, Q 2 ). 

The dedicated measurements of the low 
x physics i.e. the deep inelastic plus jet 
events, transverse energy flow in deep in- 
elastic scattering, production of jets sepa- 
rated by the large rapidity gap and dijet 
production in deep inelastic scattering are 
summarised in Fig. 4, [ flof . 

The deep inelastic lepton-hadron scat- 
tering containg a measured jet (Fig. 3a) 
can provide a very clear test of the BFKL 
dynamics at small x [|l], • The idea 
is to study deep inelastic (x,Q 2 ) events 
which contain an identified jet (aij, fey-) 
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where x << Xj and Q 2 ~ k\y Since 
we choose events with Q 2 ~ ■ the QCD 
evolution (from kj,- to Q 2 ) is neutralised 
and attention is focussed on the small x, 
or rather small x/xj behaviour. The small 
x/xj behaviour of the jet production is 
generated by the gluon radiation as shown 
in the diagram of Fig. 4a . Choosing the 
configuration Q 2 ~ k 2 -,- we eliminate by 
definition the gluon emission which corre- 
sponds to strongly ordered transverse mo- 
menta i.e. that emission which is responsi- 
ble for the QCD evolution. The measure- 
ment of jet production in this configura- 
tion may therefore test more directly the 
(x/xj)~ x behaviour which is generated by 
BFKL equation where the transerse mo- 
menta are not ordered. 

Let us now discuss the feasibility of us- 
ing the deep inelastic events which contain 
a meaured jet to identify the singular z~ x 
type of behaviour at HERA [16,41,42]. 

One practical limitation is that jets can 
only be measured if they are emitted at 
sufficiently large angles (9j > 5°) to the 
proton beam direction in the HERA labo- 
ratory frame. Large Xj jets are only emit- 
ted at small 6j ; for a given 6j we can reach 
larger Xj by observing jets with larger fcy ■ 
but with a depleted event rate. In order to 
identify the BFKL behaviour we need 
deep inelastic + jet events, with k\- ~ Q 2 , 
over an interval of z = x/xj which covers 
values of z as small as is experimentally 
possible. As a compromise we select the 
region Xj > 0.05 and x < 2 * 10 -3 . 

Fig.5 shows the predicted x dependence 
of the deep inelastic+jet cross-section rele- 
vant for HERA [16, [S]. (The cross-section 
shown in Fig.5 has been calculated assum- 
ing the following constraints: Xj > 0.05 



and Oj > 



5°, \Q 2 < k 2 Tj < 2Q 2 ). The 
continuous curves give the values of the 
cross-section when the BFKL effects are 
included. These are to to be contrasted 
with the dashed curves which show the 
values when the BFKL effect is neglected, 
that is when just the quark box approxi- 
mation is used to evaluate the correspond- 
ing differential tructure functions. The 
steep rise of the continuous curves with 
decreasing x (i.e. decreasing z = x/xj) 
reflects the z~ x BFKL effect generated by 
the gluon radiation. 

The recent HI results concerning the 
deep inelastic plus jet events which were 



reported at this Workshop [37] are consis- 
tent with the increase of the cross-section 
with decreasing x as predicted by the 
BFKL dynamics. 

Conceptually similar process is that of 
the two-jet production separated by a 
large rapidity gap Ay in hadronic colli- 
sions f43| , 44] or in photoproduction Gf| 
as illustrated in the Fig. 4c. Besides the 
characteristic exp(XAy) dependence of the 
two-jet cros-section one expects significant 
weakening of the azimuthal back-to-back 
correlations of the two jets. This is the 
direct conseqence of the absence of trans- 
verse momentum ordering along the glu- 
onic chain in the diagram of Fig. 4c. The 
experimental data on the dijet production 
have recently been obtained at the Teva- 
tron and were reported in this Workshop 



Another measurement which should be 
sensitive to the QCD dynamics at small 
x is that of the transverse energy flow in 



deep inelastic lepton scattering [47, 45]. 
The transverse energy flow in the central 
region away from the current jet and the 
proton remnants can be calculated from 
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the following formula [48|: 
dE T 



— I dkj\k 3l 



dy 

d 2 k p f d 2 k^ ( 3a s k 2 k 2 " 



vr k] 



F 2 (x/xj, k 2 Q 2 )f(xj,kl) 



(22) 



where the transverse momenta are defined 
in Fig. 6a. The function Ti describes the 
gluon radiation in the uper part of the di- 
agram in Fig. 6 a. It satisfies the BFKL 
equation as does the unintegrated gluon 
distribution /. The variable xa is re- 
lated to the rapidity y(cm) in the virtual 
photon-proton center-of-mass frame 



y{cm) = -In I 



(23) 



The rapidity y in the HERA frame is ap- 
proximately related to y(cm) by a simple 
boost 



y 



y(cm) + -In 



1. f4xE%* 



{ Q 2 



(24) 



since at small x the two frames are approx- 
imately collinear. The formula (^2|) is only 
valid in a region where both Xj and x/xj 
are suffiiently small for the BFKL equa- 
tion for f(xj, hp) and for ^(x/xj, k 2 Q 2 ) 
to be valid (we assume Xj < 10 and 
x/xj < 10 _1 ). In particular, for large Xj, 
(xj > 1CP 2 ) the more appropriate way to 
calculate the energy flow is by the use of 
the following formula: 



dy 



1 

F~2 



dk 2 u 



3a, 



7T 



^(x/xj, k 2 , Q 2 



£/.(* 



3 ) r"j , 



(25) 



which follows from the strong ordering, 
k 2 >> k 2 , at the gluon emission vertex 
(cf. Fig. 6b) where the "effective" par- 
ton combination J2 a fa can be set equal 
to g+ §(g + g)- 

The integral in the equation ( p2|) is 
weighted towards large values of the tran- 
verse momenta and so it is sensitive to the 
behaviour of the functions / and T for 
large k 2 and k 2 respectively. In the cen- 
tral rapidity region away from the current 
jet and the proton remnants the BFKL 
dynamics does predict substantial amount 
of transverse energy which grows with de- 
creasing x. 

For fixed QCD coupling a s one gets the 
x~ e increase of transverse energy with de- 
creasing x with e = (3a s /7r)2Zn2 This 
increase is closely related to the diffusion 
pattern of the solutions of the BFKL equa- 
tion for the functions / and T (cf. eq. 
@). It should be noted that without the 
Gaussian factors which broaden with de- 
creasing Xj or x/xj the integral defining 
the energy flow would be divergent. 

Numerical estimate of the energy flow 
based on the BFKL equation with running 
coupling effects incorporated confirms the 
increase of the avarage Et with decreas- 
ing x although this increase is weaker than 
for the fixed coupling case. The average 
value of Et per unit rapidity which is gen- 
erated by the BFKL dynamics has been 
estimated to be around 2 GeV. It is in- 
teresting to notice that the HI collabora- 
tion has seen the excess of transverse en- 
ergy in the forward region |37], 50 1 in com- 
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parison to the expectations based on the 
standard Monte Carlo models incorporat- 
ing conventional QCD cascades. In Fig. 7 
we confront the theoretical predictions for 
the transverse energy flow based on the 



BFKL dynamics [49] with the experimen- 
tal data. The histogram in Fig. 7 is the 
LEPTO Monte Carlo estimate of the ef- 
fects of radiation from the current jet and 
of hadronization. 

Finally let us consider the production 
of dijets close to the photon fragmenta- 
tion region and with small rapidity gap 
between the jets (Fig.4d). Absence of 
transverse momentum ordering within the 
BFKL gluon ladder modifies the theoret- 
ical expectations concerning production 
of such dijets in deep inelastic scatter- 
ing when the dominant subprocess is that 
of the virtual photon- gluon fusion. In 
particular one gets significant broadening 
of the azimuthal distributions of the two 
jets which increases with decreasing x as 
shown in Fig. 8 [51]. The distributions 



plotted in Fig. 8 are normalized to a com- 
mon maximum. The weakening of the 
back-to-back correlations between the jets 
in photoproduction has recently been dis- 
cussed in the ref. |]52[ . 

Other processes which are sensitive to 
the small x physics are the deep inelastic 
diffraction and heavy quark production. 

To sum up we have recalled in this talk 
the basic QCD expectations concerning 
the small x physics and discussed their 
possible implications for deep inelastic lep- 
ton scattering. We have limited our- 
selves to large Q 2 region where pertur- 
bative QCD is expected to be applicable. 
Specific problems of the low Q 2 , low x re- 



gion are discussed in the ref. [53]. 

We have shown that the gluon and sea 



quark distributions and so the deep inelas- 
tic structure functions should grow rapidly 
with decreasing x. This rapid type of 
increase with A ~ 1/2 is the reflection of 
the BFKL Pomeron which originates from 
the gluon radiation. 

This indefinite growth of parton distri- 
butions cannot go on forever and has to 
be eventually stopped by parton screen- 
ing which leads to the parton satura- 
tion. Most probably however the satura- 
tion limit is still irrelevant for the small 
x region which is now being probed at 
HERA. 

The observed increase of F%(x, Q 2 ) with 
decreasing x observed at HERA is in 
agreement with the predictions based on 
the BFKL equation. The explanations of 
this effect within the Altarelli-Parisi evo- 
lution equations are also possible. 

Finally we have emphasised the role of 
studying the final states in deep inelastic 
scattering at small x and have discussed 
jet production, transverse energy flow and 
dijet production. The small x dynamics 
should also show up at other semihard pro- 
cesses like heavy quark production or deep 
inelastic diffraction. 
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Fig. 1 The regions in the (Inl/x, ln(Q 2 /A 2 )) 
plane where various dynamical effects and ap- 
proximations play dominant role. 



Fig. 2 Diagrammatic representation of the 
gluon ladder contribution to deep inelastic 
scattering and of the kr factorisation formula 
©■ 
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Fig. 3 The measurements of the F2 at HERA shown together with the BFKL description 
H (continuous curves) and the MRS (A) parton analysis fit Q (dashed curves). The curves 
marked as GRV at the plots for Q 2 — 15 GeV 2 and for Q 2 — 30 GeV 2 correspond to the 
predictions of the dynamical model of ref. [p5[ . 
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Fig. 4 Diagrammatic representation of the processes testing the BFKL dynamics, (a) Deep 
inelastic scattering with the forward jet . (b) Et now in deep inelastic scattering, (c) Pro- 
duction of jets separated by the large rapidity gap Ay. (d) Dijet production in deep inelastic 
scattering. 



Fig. 5 The cross section < a >, in pb, for deep inelastic plus jet events shown as the function 
of x for three different bins of Q 2 [[l6| |42| . The continuous curves show < a > calculated with 
the inclusion of the BFKL soft gluon summation. The corresponding < a > values which were 
calculated neglecting the BFKL effects are shown as the dashed curves. 
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Fig. 6 (a) Diagrammatic representation of the formula ( p2| ) showing the gluon ladders which 
are resummed by the BFKL equation for / and T<z- (b) The representation of formula (p5|), 
which is obtained by the simplification of (|22|) when Xj is large. 



Fig. 7 The data show the Et flow accompa- 
nying deep-inelastic events with x < 10 -3 ob- 
served by the HI collaboration in the central 
region [[57J ^0|. The continuous curves show 
the BFKL predictions |||. 



Fig. 8 The azimuthal distributions of the 
jets for deep inelastic dijet events |5l| . 
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